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Abstract 

We study various matrix models with a charge-charge interaction as toy models of 
the gauge dual of the AdS black hole. These models show a continuous spectrum and 
power-law decay of correlators at late time and infinite N, implying information loss 
in this limit. At finite A^, the spectrum is discrete and correlators have recurrences, so 
there is no information loss. We study these models by a variety of techniques, such 
as Feynman graph expansion, loop equations, and sum over Young tableaux, and we 
obtain explicitly the leading corrections for the spectrum and correlators. These 

techniques are suggestive of possible dual bulk descriptions. At fixed order in the 
spectrum remains continuous and no recurrence occurs, so information loss persists. 
However, the interchange of the long-time and large- A limits is subtle and requires 
further study. 
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1 Introduction 



The black hole information paradox [T] has been a fruitful thought experiment, leading 
in particular to the discovery of gauge/gravity duality [2]. This duality in turn provides 
a nonperturbative definition of string theory as quantum gravity in AdS backgrounds in 
terms of unitary quantum mechanics, and implies that the information escapes with the 
Hawking radiation, but important questions remain. In particular, how does the argument 
for information loss, based on the low energy effective field theory of gravity, break down? 

Ref. [3], building on ideas of Refs. [11[5], presented a simple model in which this might be 
studied further. The quantum mechanical system of a single large- iV matrix oscillator and a 
single fundamental oscillator displays the key property [6] that information is lost at infinite 

but not for N finite. Since is proportional to Gn in gravity, this result demonstrates 
that quantum gravity effects are crucial to avoid information loss. In the planar limit the 
Schwinger-Dyson equation for this model closes, and moreover can be reduced to a recursion 
relation with respect to frequency [3]; analytic arguments, and numerical solution, then 
confirm the desired properties. However, a complete analytic solution seems difficult even in 
the planar limit, and a systematic study of the corrections is even more difficult. 

In the present paper we present more tractable models that have the same degrees of 
freedom, but where the previous trilinear interaction [3] is replaced by one that is quartic in 
the oscillators but quadratic in the U{N) charges. In Sec. 2 we motivate this by showing that 
the weak-coupling limit of the trilinear model leads to an effective charge-charge interaction. 
We then consider more general charge-charge interactions, and present the solution in the 
planar limit. Like the trilinear model, this displays the essential feature of a continuous spec- 
trum at infinite A^. The Schwinger-Dyson equation is algebraic, and the decay of the planar 
two-point function is power-law rather than exponential at late times. Still, information is 
lost at large A^, allowing us to address the paradox in a simpler setting. 

Our ultimate goal is to see whether the preservation of information might be reflected 
in the large-time and/or large order (in l/N"^) behaviors of the perturbative expansion. We 
are also interested in developing the analog of a bulk string/gravity description, to see how 
information is preserved in this language. To these ends we attempt to solve the model in 
several different frameworks. 

In Sec. 3 we further develop the graphical approach, beyond the planar limit. We find ex- 
plicitly the first non-planar amplitude, summing all Feynman graphs that have the topology 
of a disk with a handle. Our explicit computation of the first 1/A^^ correction demonstrates 
that the spectrum remains continuous. No recurrence occurs and information loss persists. 

In Sec. 4 we analyze the theory in terms of loop equations, recursive relations for expec- 
tation values of operators. This allows us to recover the planar term and first non-planar 
correction, and generalizes efficiently to higher orders. 



1 



In Sec. 5 we show that the correlator for this theory can be written for any in terms 
of a sum over Young tableaux. The planar limit is recovered as a large- saddle point. 

In Sec. 6 we study various models in which the matrix oscillator is generalized to a 
rectangular N x K matrix. This allows for several limits of large and small N and K where 
a more explicit solution is possible. 

In Sec. 7 we discuss the results and future directions. In particular we examine the 
possibility that one might see signs of the long-time recurrences in the behavior of the 
expansion. 

2 Charge-charge models 
2.1 Models 

First we recall the trilinear model [3]. In terms of lowering operators Aij and for the 
adjoint and fundamental, the Hamiltonian is 

H = mTr(AU) + Ma^ a + ga^ Xa + ca^ a{a^ a - 1) . (2.1) 

where X = {A + A^) I\j2m. The final term is needed to stabilize the system, but has been 
arranged to vanish in the relevant sectors a^a = 0, 1. We take M to be large so that a^a is 
essentially zero in the thermal ensemble. 
A typical Feynman graph for 

^^Mit-t') ^Ta,(t)a](t'))^ = S^,GiT,t- t') (2.2) 

is shown in Fig. [H The free thermal propagators are 

G,{T,uj) = G,{uj) = ^— , 

uo + le 

Ko{T,uj) = (— \ . 

1 — y yuj"^ — -|- te 

The only singularities in the lower-half 0^1^2,3 planes are from the adjoint propagators, at 
±m — ie, so we can evaluate these integrals by residues, leaving the fundamental propagators 
at momenta u + km for various integers k. Thus, an individual Feynman graph gives only 
poles on the real axis, as in the discussion in Ref. [5]. However, it is shown in Ref. [3] that 
this cannot be a property of the full planar propagator for any nonzero g. The point is that 
perturbation theory is singular, because higher orders of perturbation theory give higher 
order poles. The most singular graphs as u goes to its on-shell value are those in which k 
alternates between ±1 and 0, since this gives the maximum number of poles. 



) , y = e-™/^ . (2.3) 
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Figure 1: a) Typical Feynman graph for the fundamental correlator (12.21) in the trilinear 
model fl2.ip . The only poles in the lower half ti;i^2.3-planes are from the adjoint propagators. 
After evaluating residues this leaves b), where kr+i = kr ± 1. The most singular graphs 
alternate between k = and k = ±1. 



Thus we can capture the most singular graphs by integrating out the propagators at 
k = ±1, leaving an effective quartic interaction with two adjoints and two fundament alsil] 
In other words, we must do degenerate perturbation theory, because states with the same 
total number of excitations are degenerate in the free theory. Thus 



H. 



efr 



Hint ~^ 7T' Hint 

tjQ — lif) 

2im? 
_l_ 
2m? 



{a\Ai^A^,a, - a\A,^A],a,) + 0{AA, A^ A^) . 



(2.4) 



In the last line we have projected down to the a^a = 1 sector (which is annhilated by ajai), 
as relevant for G. 

The final interaction is simply a coupling of the U{N) charges of the fundamental and 
the adjoint. 

Hint = -^quQu (2.5) 
Qu = , Q^l = MjA^i - AAi ■ (2-6) 



with 



^ This is equivalent to dropping smaller terms in the recursion relations for G, as was done in [3]. 
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Figure 2: Planar Feynman graphs for the correlator (12. 2p in the model (12. 7p . The shaded 
rectangles mark the full planar propagators. Arrows point from creation operators toward 
annihilation operators. The graphs for = 0, 1, 2 are shown. 



Thus, the energy can be expressed in terms of a difference of quadratic Casimirs. This 
implies a large degeneracy; nevertheless, the model will still have enough mixing to produce 
a continuous spectrum in the large- iV limit. 

We could generalize by giving independent coefficients to the two terms in Qu. In fact, 
these two terms separately generate commuting U{N)^s, the first acting on the left index of 
A\j and the second on the right index. We obtain a slightly simpler model (in terms of its 
Schwinger-Dyson equation) by keeping only one term, 

Hir^t = -hqi.Qu , Qu = A] A . (2.7) 



2.2 Planar solution 



Now let us solve these in the planar approximation. For the qQ model (12.71) the nontrivial 
planar graphs involve a cycle with n + 1 vertices, giving the Schwinger-Dyson equation shown 
in Fig. El Thus, 



n=0 



Lo(r, -ui) J]^ \ G{T,uj - uji+i - ui)Lo(T,uJi+i) ,(2.8) 



(27r)" 



1=1 
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where Lq[T,uj) is the thermal AA^ propagator 

Lo{T, uj) = -^ ( . (2.9) 

At large M, G{T,u!) has singularities only in the lower half tu-plane, and Go and Lq both 
fall as l/u at large frequency, so we can close the Ui integrals in the lower half-plane and 
pick up residues only from Lq. This gives 

S4T,u;) = y{^^] GiT^uf . (2.10) 



1-2/ 



The Schwinger-Dyson equation becomes 



- iujG{T, u;) = l-- ~ , (2.11) 

I — y + iXG[T, uj) 

or 

- uXG^{T, Lu) + i{l -y){uj + \)G{T, u) + {1 - y) = , (2.12) 
where A = hN is the 't Hooft coupling. The solution is 

G{T,u) = %^(x + u;-y/{uj-uj+){uj-uj^)) , = X^-Ll^Nl , (2.13) 
2ujA \ / 1 — 2/ 

This has a pole of spectral weight 1 — y at u = 0, and a cut from U- to 0;+. 

For the original model fl2.5p the only modification is the inclusion of an additional copy 
of each cycle but with the arrows reversed, so that 

n+l 

with A = hN. Then 



S„(r,^) = !!,-(-!,)''] I -^1 G(T,u,r. (2.14) 



l-y + i\GiT,uj) l-y-iy\G{T,uo) 

This reproduces the cubic equation for G(T, cu), Eq. (32) of [3j, which was obtained from the 
weak-coupling approximation to the recursion relation of the trilinear model. Again, there 
is a branch cut on a finite segment of the real axis. 

2.3 Black hole physics 

The continuous spectrum implied by the cut in these models is the signature of a horizon. 
As noted in Ref. [3], the cut is absent at zero temperature, and also below the Hawking- 
Page transition (which we can simulate in this one-matrix model by imposing the singlet 



5 



constraint). The Fourier transform of the cut gives a behavior at late times. Although 
this falls off more slowly than the exponential for the real black hole, it is still inconsistent 
with the properties of a system with finite entropy, and so with the exact correlator. Indeed, 
the energies in this model are all multiples of h = X/N, so there are regular recurrences with 
period 2nN/X. This can be seen by writing the interaction in terms of quadratic Casimirs, 
as we will do in Sec. 5.1. This time is much shorter than for a fully thermalized system, 
for which the energy splittings are of order e~^^, but still presents us with a version of the 
information paradox. We should note that the more general model 



cannot be written in terms of commuting Casimirs for generic hi ^2, and so may give a more 
realistic model of the black hole. 

3 Non-planar corrections 
3.1 Schwinger-Dyson equation 

We now consider the full Schwinger-Dyson equation, including non-planar corrections. We 
focus henceforth on the qQ model (12.71) . It is useful first to carry out all of the loop integra- 
tions, as in the previous discussion. The number of loops is equal to the number of adjoint 
propagators, so we can take the adjoint propagator momenta Ui as integration variables. 
We orient these in the direction of the arrow on the fundamental propagator, as in Figs. 1, 
2. The momentum on the fundamental propagator therefore always involves —Wi, and so 
this propagator contributes no poles in the lower-half Ui plane. Thus we can close the loop 
integrals in this half-plane picking up the pole at m — ie for each forward propagator (one 
whose arrow is parallel to that on the fundamental line) and at — m — ie for each backward 
propagator (antiparallel to the fundamental line). Further, since each vertex contains one A 
and one A\ there are always equal numbers of +m and —m loop momenta flowing on any 
internal fundamental line, so this is always at Ci;intcrnai = i^extcmai- This is simply a repetition 
of the point made in Sec. 2.1, that this model isolates the propagators with k = 0. 
We therefore have the Feynman rules 



int 



(2.16) 



^ '^external ; 



1 




(3.1) 
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The evaluation of the amphtudes is thus reduced to counting graphs, 

graphs 

where v is the number of vertices, b is the number of backward propagators, and g is the 
genus of the graph. 

We now reduce the sum over graphs by summing over certain classes of subgraph. 
Start with the sum of the one-particle- irreducible (IPI) fundamental self-energy graphs 
S{T, u, h, Go), in terms of which the full Schwinger- Dyson equation reads 

G{T, u) = Goico) + Go{uj)S{T, cu, h, Go)G{T, cu). (3.3) 

We denote explicitly its dependence on the coupling and the functional form of the propagator 
for the fundamentals. Consider self-energy corrections on internal propagators. Because we 
are assuming that the fundamental oscillator mass M is large, there are no loops of the 
fundamental field, and so there are no corrections to the adjoint propagator, only to the 
fundamental propagator. By packaging corrections on internal fundamental propagator as 
G, we only have to sum over graphs without them. What we get is precisely the two-particle- 
irreducible (2PI) fundamental self-energy S'2Pil^ so 

S{T,u,h,Go) = S2MT,uJ,h,G). (3.4) 

Next, for any vertex, there is a geometric series of graphs obtained by expanding it into 
n + 1 vertices, each connected to the next by a forward adjoint propagator and a fundamental 
propagator. This is illustrated in the SD equation of Fig. 2, where all of the interaction 
terms can be obtained from the n = term in this way. Thus, we can sum all such series 
into effective vertices. We will refer to the propagators to be summed as trivial forward 
propagators, a term which we will explain in Sec. 13. 2[ So we can restrict the sum to graphs 
with no trivial forward propagators (NTF) but with the geometric series incorporated into 
the vertices: 

S{T,u;,h,Go) = Sr4T,u,h',G) , h' = -. (3.5) 

Here Si* denotes the sum of the graphs that are fully irreducible (I*), i.e. IPI, 2PI, and NTF. 



2PI graph is one that cannot be separated into two pieces by cutting two propagators. Such a graph 
clearly cannot have an internal propagator correction. If a graph in S can be separated in this way, at least 
one fundamental propagator has to be cut, because all vertices lie along a single fundamental line. Since 
the number of forward minus backward propagators emerging from any subgraph is always zero, the second 
propagator to be cut must be a fundamental line, isolating the propagator corrections in between. 
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Figure 3: Schwinger-Dyson equation for the full propagator. The sum runs over all I* graphs. 
Lines with shaded semicircles denote full propagators. Vertices with shaded semicircles 
contain the summed trivial forward propagators. When we refer to the genus of this graph 
(as in Sec. 3.2, where we sum all genus one graphs) we mean the explicit genus, not taking 
into account the non-planar corrections implicit in the shaded circles. 

One form for the Schwinger-Dyson equation is then 

G(T, u) = - + -Si, (T, u, h', G) G(T, u) . (3.6) 

The total contribution of a given I* graph to Si* (T, tu, /i', G) G{T, u) is 

^y^V-^s , w = -ih'NG{T, uj)/{l - y) . (3.7) 

The first few terms are shown in Fig. 3. The planar contributions have been collapsed into a 
single term by the summation over graphs containing trivial forward propagators, but even 
at the first non-planar order, the number of I* graphs is infinite. Summing these is our next 
exercise. 

3.2 The l/iV^ correction 

If we take the trace of the fundamental propagator, which is just NG(T, u), then the planar 
graphs are those that can be drawn on a disk, while the A^~^^ corrections come from graphs 
that can be drawn on a disk with g handles [7]. Thus we wish to enumerate all fully reducible 
I* graphs that can be drawn on a disk with one handle (Fig. 4). We have marked with an 
X the point where the ends of the fundamental propagator are joined, because this enters 
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Figure 4: Disk with one handle. The x indicates the point where the ends of the fundamental 
line are joined. The dashed lines mark the A and B cycles. 

into the Feynman rules. We have also marked A and B cycles; every adjoint propagator 
is homotopic to pA + qB for some integers p and q. For the present discussion we are not 
distinguishing an orientation on the adjoint propagators. 

The possible trivial propagators, p = q = 0, are very limited in a 2PI graph. As is 
clear from Fig. 2, the forward propagators that have been summed into the NTF vertices 
are trivial, and in fact these are the only trivial forward propagators as we have already 
mentioned. To see this, note that a contractible forward propagator divides the Riemann 
surface into two pieces (that which it crosses when contracting, and that which it does not). 
It follows that if we cut the fundamental propagators that attach to each end of this forward 
propagator, the surface is separated into two. Because the graph is 2PI, this is only possible 
if the two ends are connected by bare fundamental propagator. Thus it is of precisely the 
type summed by the NTF condition. 

If a backward propagator is trivial, a similar argument shows that it also divides the 
Riemann surface, and so would cutting on the adjacent fundamental lines. This is also 
excluded by the 2PI condition, the only exception being when the ends of the backward 
propagator are the first and last adjoint propagators to attach to the fundamental line, since 
the IPI condition omits these adjacent propagators. For example, all the planar graphs in 
Fig. 2 have such a trivial backward propagator. So for any non-planar graph, either all 
propagators are homotopically nontrivial, or there is a single trivial backward propagator 
which separates the x from the rest of the graph. 

To enumerate the nontrivial propagators, we must be careful not to overcount, because 
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the SL{2, Z) modular group of the torus allows us to draw the same graph in different 
ways. Therefore, we specify that as we move along the fundamental line in the direction of 
the arrow, the first nontrivial propagator that we encounter is homotopic to the A cyclejf] 
Similarly, we specify that the second nontrivial propagator that we meet (excluding those 
homotopic to the first) is homotopic to the B cycle. This fixes the modular group, and so 
we can count freely. Note that there is at least one propagator along the A cycle and one 
along the B cycle, or else the graph is actually planar. 
There are five possible kinds of propagator: 

1. Ui propagators which go under the handle from the left of the marked point {A cycle 
propagators). As noted above, ni > 1. 

2. 77-2 propagators which go directly along the handle (the B cycle). Again, ri2 > 1 if the 
graph is non-planar. 

3. propagators which go along the handle in a "twisted" way. These are homotopic to 
A + B. There are genus-one graphs without these twisted propagators, so > 0. 

4. n4 propagators which go under the handle from the right of the marked point. Again, 
these are homotopic to the A cycle, but there need not be propagators of this in this 
group: ^4 > 0. 

5. trivial backward propagators as discussed above, where G {0, 1}. 

In Fig. 5 we depict in various ways the case rii = 2, n2 = = = = 1. 

The integers rii are subject to several parity constraints. Because the operators A and A^ 
alternate along the fundamental line, there must be an even number of propagators attaching 
between the ends of any given propagator. By applying this to given propagators of types 1, 
and 2, we find that n2 + n-s and ni + n3 + 724 must be even. Type 3 propagators give the sum 
of these conditions, type 4 give the same condition as type 1, and type 5 give no conditions. 
This leaves four cases: 

(ni,n2,n3,n4) = OOOE, EOOO, OEEO, EEEE , (3.8) 

{E = even, O = odd), and in each case ns = or 1. 

It remains to determine the number of backward propagators b, since each of these brings 
in a factor of y. The number is approximately half the total number of propagators, but a 
precise count requires us to examine separately each parity case (13.81) . For ns = one finds 
respectively 

6=(n + l)/2, (n + l)/2, (n + 2)/2 , n/2 , (3.9) 

■^We can exclude a integer multiple of A (which would not be modular-equivalent to A) because such a 
propagator would intersect itself. 
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Figure 5: a) The I* graph rii = 2, n2 = = n4 = = 1 drawn on the disk with handle, 
b) The same graph drawn in the usual double line notation, c) The same graph drawn as a 
torus with a hole. 
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where n = ni + n2 + ns + 77,4. For 725 = 1, 



b = {n + l)/2, (n + l)/2, n/2 , (n + 2)/2 . (3.10) 

Recall that the weight (13. 7p of a given graph is w^y'', with v = n + n^. We sum over the 
integers (rii, 71.2, ris, 724) with the limits rii > 1, ^2 > 1, > 0, 77,4 > 0, and also over n^, 
separating according to parity case: 

OOOE EOOO 



tf^|/^(l + tf )^(1 + Wl/) 

(1 - w^yy 
x^y^(l — a;)^(l — x[l — y]) 



Here x = w/(l + w) = -i\G{T, to) /{I - y) = -i\yG{T, cu). 
Expand the correlator in l/N"^, 

GiT,u;) = G^'\T,u) + l-^G^'\T,u;) + o(^-^^ , (3.12) 

where G'^°^(T, u;) is identified with fl2.13p . Also we define Xq = —iXyG^^\uj). The l/N"^ term 
in the SD equation is then 

= 1 (-a,G(i)9,[5i.G],=o + [SiA9=i) , (3.13) 
where [Si*G]g=Q = wy = xy/{l — x) from the single planar I* graph. The solution is 

rW/T ,] = iy^xljl - xo)\l - xo[l - y]) 

^' ^ {1 - 2X0 + xl[l - y]nu[l - xof - Xyy) ■ ^' ' 

The RHS of Eq. (13.141) is a rational function of u and of xq, so its branch cut comes 
only from that of Xq: it is in the same place as for the planar amplitude G^^\T, lu). For real 
uj, IieG^^\T,uj) is non-zero only if ReG*^°-'(T, u;) is also non-zero. Furthermore the contin- 
uous spectrum of Re G'-''-* (T, cj) is not modified by the leading perturbative correction 
Re G*^^^(T, ci;). Note that from (13.61) . (13.71) . each I* graph contributing to Si*G has the same 
branch points as the planar amplitude, and the same continuous spectrum. This need not 
hold after summing an infinite series of I * graphs, but we have found that [Si*G]g=i (13.111) 
is a rational function of Xq and uj, and so does not introduce new cuts. We expect that this 
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will continue to hold at higher genus, so the continuity of the spectrum and the positions of 
branch points will be the same at any finite order in 

However, even this first correction does change the nature of the branch point. By 
inserting the result (12.131) for the planar propagator, one finds that 1 — 2xo + Xq[1 — vanishes 
as (a; — uj±Y^'^ at the ends of the cut, so that the denominator vanishes as {to — uj±Y. The 
branch point behavior of G^^^ is then more singular than that of G^'^^: there is an (cu — uj±y 
double pole, as well as a subleading {to — uj±)~^^'^ branch cut, as opposed to (a; — co±y^'^ in 
the planar term. (This effect, which seems rather accidental in the present approach, will be 
more evident in the next section.) This implies that G^^^/G^^^ grows at long times. In the 
conclusions we will discuss the possible relevance of this for the information problem. 

It would be very useful to extend the graphical solution to all orders in l/N"^. The main 
challenge seems to be the treatment of the modular group. It would be interesting to find 
in particular some string field description, in which the sums over adjoint propagators can 
be thought of as a string propagator. We hope to return to this in future work. 

4 Loop approach 

One of our goals is to find a 'bulk' description of our model, one analogous to the gravity 
side of gauge/gravity duality. Bulk quantities are invariant under the gauge symmetries 
of the CFT, so we must work with U{N) invariant objects. Recasting gauge theory in 
terms of invariants has long been a tantalizing idea for connecting gauge theory with string 
theory [8l[9], but it has been difficult to implement. In our model, it turns out to be a useful 
way to calculate. 

The invariants that we work with are TrQ", and generating functions for these. We have 
to be careful about ordering because Qij = A^^^A^j is a matrix whose elements are operators. 
We specify the natural ordering whenever there is a matrix multiplication, e.g. 

Qlj = QikQkl ■ ■ ■ Qmj ■ (4.1) 

That is, we think of Q as acting on the tensor product of the index space and the adjoint 
Hilbert space. We will encounter expressions where other orderings arise, and for these we 
write the indices explicitly. Also we define Qij = AikA\-. Note that the matrix elements 
commute with the Qki- 

Because the thermal ensemble has no fundamental excitation, we can write the correlator 
(O as 

NG{t) = e{t) ^aie-^'^"l^^'="^*aj^^ . (4.2) 

(We will write {■ ■ ■)t simply as (...) hereafter.) Since we consider the limit that the mass 
of fundamental field is very large, the number of fundamental field is always one. If we 
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expand the exponential, only adjacent oscillators can be contracted. Thus we can express 
the correlator in terms of invariants as 



NG{t) = 9{t) {Ti e~''"^') , 

The trace is over the A^-dimensional matrix indices, while the (thermal) expectation value 
is in the free adjoint Hilbert space. 

4.1 Thermal loop equations 

In a thermal ensemble, the Heisenberg operators satisfy 

{A,i-tP)O,,i0)) = {0,mAjm ' (4-4) 
for any Oji. The free field equation Aij{t) = e~*"**yljj(0) then implies 

{0,,A,) = y {A,0,,) , (4.5) 

or 

(0,,A,,) = ^([A,„0,,]) (4.6) 

We now omit the common time argument from all operators. The 'loop equation' (14.61) de- 
termines all expectation values iteratively, since the right side has two fewer mode operators 
than the left. 

Let us first illustrate this with some simple examples. For Oji = Aj- we obtain from 
Eq. gl]) 

{TtQ) = -^. (4.7) 
1 



For Oj, = {A^AA^),i = (QA^), 

• (d^O]\ = /Tr (9.0) /V2\ = ^'^y 

y) 



(Tr Q^> - ^ (Tr (g ^- Q)) ^ ^ (Tr (2g) + N^) ^ . (4. 



For O,, = QkkA], 

1-y \ I l-y iX-y) 



(TrQTrQ) = (tiQ) + ^ (TrQ) = + (TrQ)^ . (4.9) 



In particular we obtain the connected contribution, which is down by 0{N 



2\ 
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Now let Oji = {Q^A^)ji; this will require a bit more work. First, the loop equation gives 

n 

(TrQ"+^) = ^^^{^TrgPTrg'^-P^ . (4.10) 



If Q were a c-number matrix, the cyclic property of the trace would give Tr = TtQ^, but 
here we have commutators]^ 



r-l 

Tr g'- = AjQ'-'Al = Tr g^ + ^ Tr Q" Tr g^"""^ . (4.12) 

g=0 

To solve this, we introduce the resolvents 

/(^) = 7T^, f{z) = J—; F{z)=TTf{z), F{z) = TTf{z). (4.13) 
^ - z Q - z 

Multiplying Eq. (14.121) by —z~^~^ and summing r from to oo gives 

F{z) = F{z) - F{z)F{z) , (4.14) 

so that 

F(^) - ^ . (4.15) 
To apply this to the recursion (14.10p . multiply by —z'"^'"^ and sum, giving 

w + m - -y^ {Fi^m^)) - (y^) . (4.16) 

This can be used immediately to obtain the planar propagator, but we first derive further 
results that will be useful beyond planar order. 

The next case is Oji = (TtQ"') {QPA^)ji. The loop equation is 

^(Trg'^TrgP+i) (Trg^Trg^Tr gP-^) + ([A,„Tr g"](gP At),,) . (4.17) 



Alternatively, one can compute (TrQ"'''^); and hence the correlator (|4.3p . by applying both (|4.6p and 

(Or^4) = T^([^^-4]) (4-11) 

recursively to eliminate the rightmost oscillator at each step. After all oscillators are contracted, the re- 
sulting terms are in one-to-one correspondence with the Feynman graphs for G. Equations (|4.6p and (|4.1ip 
respectively generate backward and forward propagators, and reproduce the Feynman rules The ends 

of propagators in a Feynman graph appear in the opposite order to how the corresponding oscillators are 
inserted. 
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To simplify the last term we use 

n-l 

lik 



q=0 
n-l 

q=0 

= E Q^r^'QlA^ ■ (4-18) 

g=0 

We have used the fact that AQ = AA^A = QA, and that the matrix elements of Q and Q 
commute. If we could reverse the order of Qlj and Aik this would simplify. Thus we proceed 

q-l 

Ql.A, = {AQ% - ^(TrgO(^Q'"^"')o- • (4-19) 

r=0 

Now multiply by —z'^'"^ and sum. Then 

f{z)kjAk = (1 + F{z)){Af{z))ij . (4.20) 

To use this, multiply Eq. fl4.18p by — 2^^"" and sum on n, 

[Aj,F{z)] = -f{z)uf{z)kjAik 

= -~f{zUl + F{z)){Af{z))i, 

= -(l + F(^))(/»M),, , (4.21) 

so 

[A,,F{z)]{Q^A% = -(1 + F(.))Tr(/(z)2Q^+i) . (4.22) 
Now multiply the relation fl4.17p by and sum from {n,p) = (0,0) to obtain 

^ {F{z){F{w) + N/w)) = - (^F{z)F{w)F{w)) - ((1 + F(z))Tr (f (z)[/» + 1/w]) 

(4.23) 

We can simplify further by using partial fractions, 

Pm'M + IH = 9. ( + M I . ,4.24) 

\ z — w w 

Then 

^—^{F(z)(wF(w) + N)) = -w (f(z)F(w)F(w))-( (1 + F(z))d, I ^-^(^) " ^-^(^ ) 
y \ / \ \ z — w 

(4.25) 
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Finally, we consider a very general case Oji = F{zi)F{z2) ■ ■ ■ F{zi){f{w)A^)ji. The com- 
mutator ( ]4.2ip can be written 

A,F{z) = {J'{z)A),, , J^{z)ki = F{z)6ki - (1 + F{z))p{z)ki • (4.26) 

Then {OjiAij) = y {AijOji) becomes 

{F{z,)F{z2) . . . F{zi){wF{w) + N))=y (^A,F{z^)F{z2) . . . F{zi){A^ f{w)),,) 

= y (Tr^(zi)^(^2) . ..Hzi){wf{w) + 1)) . (4.27) 

This is our most general form for the loop equation. The trace acts on the matrix indices 
of JF and /. Note that the matrix elements of f{zi) commute with F{zj) and so can all be 
brought to the right. 

By partial fractions the RHS can again be written in terms of expectation vales of F, 
and their derivatives. For example, in the case that zi = Z2 = ■ ■ ■ = zi, 

(FizYiwFiw) + N))=yJ2(l) (-1)"^ (f(^)'-'"(1 + F{z)r Tr fizf^^iwRw) + 1)) 

I 

= y(^F{zy{wF{w) + N)'j+yY, 

m=l 

(F{zy~-{1 + F{z)rdl-'' |^ ^%)-^F(^) j \^ ^428) 

Here we used 

p-{z) [wfH + 1] = T^^^^r-^ f zm-wf{w) \ 

[2m — 1)! \ z — w I 

4.2 expansion 

We have not yet made any use of large A^. In the 't Hooft limit, Q,z,w are of order A^, 
F is of order 1, and / is of order A^~^. We therefore rewrite the loop equation in terms of 
0(m) = Nf{Nu), $(m) = F{Nu) and 4>(m) = F{Nu). Note that 

J^{Nu)ki = Hu)6,i - ^(1 + ^u))4>''{u)ki , (4.30) 

so the second term is nonplanar. Therefore from Eq. fl4.27p 

mm)^U2) . . . ^ui){v^v) + l))=y ($(ni)$(n2) . . . ^ui){v^v) + 1)) + Oil/N^) , 

(4.31) 



-1^ 



m J (2m — 1)! 
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or, using $ = $/(l + $), 

. . . [v^v) + l-y- 1/^YT$m}) ° • ^^'^^^ 
In the planar limit the expectation value factorizes, and so we have 

v^liv) + (1 - y)(l + v)%{v) + (1 - y) = , (4.33) 

where ^o{v) = ('^'('y))pianar- T^^^^ reproduces the planar result fl2.12p for ^o{v) = i\G(T,v\). 

We have assumed large-A^ factorization, but we should be able to derive it from the 
loop equations, since we have argued that these are complete. We give a somewhat formal 
derivation, as follows. In Eq. (14.321) let the Ui all be equal to v. By forming a power series 
we can conclude for any analytic function r($(t;)) that 

r($(^)) i^v^v) + l-y- y^YT^^]) ° • ^^-^^^ 

The thermal ensemble produces some probability distribution for $(t')lf| Since Eq. (14.341) 
holds for arbitrary r, it must be that the distribution is concentrated on the zeros of the 
expression in the bracket. There are two such zeros, but these have different large v be- 
haviors, 0{l/v) and 0(1). Using the known asymptotic behavior we can conclude that the 
distribution is a delta function on the first solution. 

We can now solve iteratively for the higher orders, expanding the general loop equa- 
tion (lOTI) as 

$(mi)<I>(m2) . . . $(««) < v^{v) + l~y-yv 



1 + <^{v) 



i=l \ * 
l<i<7</ ^ 



{Ui - Uj){Ui - V) {Uj - Ui){Uj - v) {V-Ui){v-Uj) 

+0 ( ^) • (4.35) 



^ Since is the generating function of all Casimirs TrQ", this is precisely the probability distribution 
of Young tableaux discussed in Sec. [Sj It is concentrated on a typical tableau determined there. The explicit 
formulas for Tr Q" in terms of tableaux can be found in [TO] . 
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To obtain the propagator to 0{1/N'^) it suffices to take all Ui to v as in Eq. fl4.34p . giving 

[^v)-%{v)] [v^v)-{l-y)/%{v)\ 



1 + ^{v 
y I dr 



{^[v)){\ + ^{v))(2^\v) + v^"[v)) ) (4.36) 



2iV2 

The LHS is the same as in Eq. (14.341) . but now written in terms of the planar solution 

$o(-) = (l + - - v/(- -.,)(.-._)) , = i±l±M , (4.37) 

2v \ / i- — y 

by using Eq. (14.331) . To solve this to order l/N"^, simply choose 

1 + $ 
- (1 - y)/^o 

Inserting the planar solution ^o{v) on the RHS, Eq. (14.361) becomes 

^(^) = Mv) - ^^(*o(^))(l + + v^'^iv)) + O (ij) . (4.39) 

We can eliminate $o in favor of $o using $o = $o/(l + '^o)- Now derivatives of the quadratic 
equation (14.331) allow us to eliminate $q and $9. The 0(1/A^2) term in (14.391) is then a 
rational function of $0; and its equality with graphical approach result (13.141) can be shown 
by using the quadratic equation again. 

This representation of the solution makes evident the fact that the branch cut in the 
propagator becomes more singular: the second derivative converts the {to — uj±Y^'^ to (a; — 
uj±)-'^/^. This still misses the full singularity, as there is an additional {u — uj±) ' hidden 
in dr/d^. It is also evident, from the form of the higher nonplanar terms, that there there 
will be two additional derivatives at each further order in and so the singularity will 

become two powers worse at each successive order. Therefore at late times, the non-planar 
corrections G^^\t) grow as t^g-s/a^ 



5 Solution by summing over tableaux 

For a space-time picture to emerge from a quantum mechanical model, an important step 
would be to find variables that are to be path-integrated and identified with the metric. 
It is conceivable that such variables are in fact discrete at the fundamental level, forming 
continuum approximately only in the large limit. Here we use group theory and propose 
one such description for the qQ model (12.71) : the sum over Young tableaux replaces the path- 
integral over the metric. In the large N limit, we recover the expression for G obtained in 
Sec. 2 and 4. 
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5.1 Sum over tableaux 

We start with the Fourier transform of (14.21) 

where q ■ Q = QijQji- By inserting the completeness relation in the middle, we can write the 
correlator as a sum 

-.GM = l(l-y)-X:(^|;^^|B). (5,2) 

The sum is over the states |-B) with one fundamental excitation {a^a = 1). Such states form 
a large representation of the U (N) x U (N) acting on the left and right indices of A^; we now 
decompose it into irreducible ones. 

The states \B) are obtained by acting one a| and several Alj, so they form the represen- 
tation 

(□,1) ® © Sym' □,□ (5.3) 

1=0 

of U{N) X U{N). We symmetrize because are bosonic. To decompose the symmetric 
product we compute the character, for (f/, V) G U{N) x U{N). Denoting by k = {ki, k2, . . .) 
an infinite series of non-negative integers, we have 

Trgyj^!(g^|-I)(f/ ^,V) = ^ -^{0\A^ni ■ ■ ■ Ai^mUnimi ■ ■ ■ UniniiVj^i^ . . . Vj^^Al^^j^ . . . ^j^^^-JO) 



I 



= J2Il-^(^'U')''(^'^')''- (5-4) 

To rewrite this, recall that the vector k labels a conjugacy class C{k) of the permutation 
group Si such that / = '^jkj. Each irreducible representation of Si is labeled by a Young 
tableau R whose number of boxes \R\ equals /. The tableau also specifies an irreducible 
representation of U{N). The Frobenius relation 

HiTrWp = J2xR{C{k))TTRU (5.5) 
j R 

and the orthogonality 
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where XR{C{k)) is a character for Si and Tr^f/ is the trace of U in the representation of 
U{N) specified by R, allow us to rewrite ( 15. 4p as 

J2^T,^rn^im(U'\V) = J2'^TjiUTTjiV = J2'^TT^U-'TTnV. (5.7) 

I R R 

Thus 

© Sym' (□,□) = ©(i?,i?). (5.8) 

/=0 R 

The sum is over all Young tableaux with at most rows. Let us recall that a Young tableau 
R represents a partition I = Ri + R2 + . . . + Rn of an integer / = \R\ by a non-decreasing 
no n- negative integers -Ri > -R2 > • • • -Rat > 0; is visualized by a collection of boxes where the 
z-th row has Ri boxes; and specifies an irreducible representation of U{N) as in (15. 8p . The 
box in the z-th row and j-th column will be denoted by 

We further decompose D® -R using the Littlewood rule. Let us consider a large rectangle 
that contains R and has rows. The complement of R is another Young tableau. The 
Littlewood rule instructs us to add a box to the complement so that the result is again a 
Young tableau. This is equivalent to removing a box [i, Ri) at an outward corner of R, 
resulting in a Young tableau R' = R — (i, Ri) with R'j = Rj — 6ji. For i = N, we relax the 
conditions on R' so that i?'^ can be —1. Then 

n0R= © {R-{t,Ri)). (5.9) 

(j,-f?i):corner 

Thus, the Hilbert space decomposes as 

(□,l)©©(^,/2) = © © {R- {i,Ri),R). (5.10) 

R (i,/?i):corncr 

In the first expression, the generators q and Q act only on □ and R, respectively. In the 
second, q + Q act on i? — (i, Ri) = R . Writing the interaction in terms of the quadratic 
Casimirs 

-hq-Q = hTT{q^ + Q^-{q + Q)^) (5.11) 
and using their values [lOl [11] 

N 

Trg^h = N, TtQ% = Rj{R, - 2j + 1 + N), 

Tr(g + Q)% = R'.iR'j " 2j + 1 + N), (5.12) 
j 
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we find tlie energy spectrum —hq ■ Q = h(N — i + Ri). Noting TrQ|;^ = \R\, we can finally 
rewrite ( 15. 2p as 

-tG{ij) = (1 -y)^'^yl^l(dimi?)2 n{uj) = {n{uj)), (5.13) 

R 

where we have defined 

^ ^ LU-h{N-t + Ri) 

('i,Ri):coincT 

^^=iV di^ ' ^« = 1- (5-15) 

(i,-Ri):corncr 

In this form, the correlator is given as the average of the "resolvent" in the statistical 

ensemble of Young tableaux with weights oc ?/l^l(dimi?)^. 

5.2 Typical tableau at large 

Statistics of tableaux with various weights have been studied in the mathematics litera- 
ture. See [12] and references therein. In physics, similar problems appeared in studying the 
Seiberg-Witten theory using instantons [13] . 

As shown in Fig. [6l each tableau R specifies a piecewise linear curve y = f{x) such that 
f{x) > \x — 1|. Since R can have at most rows, we must have f{x) = —x + 1 for x < 0. 
The number of boxes |-R| is proportional to the area of the region \x — 1\ < y < f{x). The 
dimension of the U{N) representation is given by the formula [H] 

where the hook length h{i,j) for the box (i, j) is defined as 

hit,j)=R,-j + Rj -t + 1 (5.17) 

with Ri and Rj being the number of boxes in the z-th row and j-th column, respectively. 
The box is mapped to {x,y) = (1 + (j — i)/N, (j + i)/N), so when is large, dimi? is 
approximately given by 



iV2 

dimi? = exp I dxdy{\ogx — \og{v — u)). 

'\x-l\<y<f{x) 



2 



The quantities u and v are given hj u = x — f{u) + y and v = x + f{v) — y, and are indicated 
in Fig. [61 Thus the weight in fl5.13p is now given by 
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(0,0) n I V X 



Figure 6: The Young diagram i?, shown rotated and inverted, defines a piecewise hnear curve 
y = f{x). The upper corner of the box (z, j) is placed at {x, ?/) = (! + (j — i)/N, (j + i)/N). 
The hook length is proportional to v — u. 

yl«l(dimi?)2 = e-^'^[^l, (5.19) 
where we have defined the "action" S[f] by 

2^[/] = jdx {^-2\ogx) (/(x)-|x-l|) 

+ I dudv{l + f{u)){l- f{v))\og{v ~u). (5.20) 

In the large limit, the ensemble is expected to be dominated by the tableaux whose 
shapes are almost identical. In other words, the probability distribution for / should be 
sharply peaked around some /*, which specifies the shape of a typical tableau. We assume 
that f^:{x) is convex, so that —1 < fi{x) < +1 for X- < x < x.^ and fi{x) = ±1 for x ^ x±. 
The values of x± will be determined below. 

To find as a minimum of S[f], let us take the variation 



1 /"°° /m 
^S[f] = -- J dx6f{x)(^j;X-2x\ogx + 2x + I[f]{, 

- j rfa;5/(x)(--21ogx + /[/r(x)), (5.21) 



2 

where we have defined the integral transform 



J[/](x) = / dx ((x — x) log \x — x'\ — (x — x')) f"{x). (5.22) 

J — oo 

We see that the function /^(a;) must satisfy the following conditions to be a minimum of 
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a) 



-1 +1 v^i(oo) 

b) 



-1 



c) 



Figure 7: a) The upper half 2;-plane on which the function (p{z) is defined, b) An auxihary 
half-plane we introduce to determine (p{z). It is useful to map the z-plane to this half w- 
plane by a Mobius transformation, bringing x± to canonical locations ±1. c) The image of 
the half z-plane by ip. 



S[f]: 



777 r 

--21ogx + /[/.r(x)| 



> when x > x+ or x < X-, 
= when x_ < x < x+. 



(5.23) 



For to be a minimum, 5S has to be zero or positive under the variation 5f which is 
constrained as Sf{x) > when f'{x) = ±1, leading to the inequality. 

To find such consider an analytic function (p{z) defined on the upper half-plane as 



I 

IT 



m 
T 



2 log z + 



dx log{z - x)f"{x) 



(5.24) 



Since log(x + ie) = log \x\ + 7xi6{—x)^ the restriction to the real axis is 

% / TTl \ 

ip{x) = lim ifix + te) = /:(x) + 2ei-x) + - /[/.]'(x) + - - 2 log |x| . (5.25) 

e~*0 71 \ i J 



Then we see from fl5.23p that the conformal map ip takes the interval x+] to [—1,1]. The 
interval [0,x_] is mapped to the vertical line passing the point z = —1, while the half lines 
[a;+,-|-oo) and (— oo,0] are are mapped to the vertical line passing z = 1. Because of the 



logarithm we have \imip[±e) = ±1 -|- ioo, and the large |a;| behavior 

e— >0 



/[/,]'(x) = 21og|a;|--+fl ^'^ 



1 / 1 
-^ + o — 



implies that 



m 2il /I 

ip{z) = 1 + i— ^ O { — ] as z ^ oc. 

vrT n z V ^ 



(5.26) 



(5.27) 



In particular, 93(00) = 1 + im/'n-T. 

These conditions can be visualized by saying that ip maps the upper half-plane in Fig. 
[7^.) to the semi-infinite strip in Fig. [Tt). The four points {0,x±, 00} should be mapped to 
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{cxD,±l,l + im/nT} respectively. To determine ip, it is useful to first consider a Mobius 
transformation that maps the three points and x± to oo and ±1 (see Fig. Eb)) 

^'^'^ = ix,-x.)z • ^^-2^) 

The conformal transformation from the half-plane to the strip that respectively maps oo and 
±1 to +ioo and ±1 is uniquely given by (p2{'w) = (2/7r) arcsin(ty). Therefore we conclude 
that if is given by the composition ip = ° Vi- 

We still need to determine x±. The map of the fourth point (oo) imposes one condition 
on x±: 



sm 







E(- 





y,,(oo) = ^+±^. (5.29) 



We get the other condition by comparing the expansion of (p2°<Pi around z = oo with fl5.27p . 
Equating the 0(1/2;) terms, we find that 

x+x^ = 1. (5.30) 

Thus x± are given by 

x+ = coth— -, a;_ = tanh— -, (5.31) 

and the conformal map if by 

ifiiz) = — arcsin (cosh sinh | . (5.32) 

Finally, the dominant tableau has the shape given by 

flix) = — arcsin (cosh — sinh — - | (5.33) 

•'^^ ^ vr \ 2T X 2T J ^ ^ 

within the interval [x„,x+], and f{x) = |x — 1| outside. If one uses this to compute the 
0{l/x'^) term in /[/*]'(x) and compare it with (I5.26p . one finds that |i?|/A^^ = y/{l — y) as 
expected from Bose statistics. 

It is also possible to solve the conditions fl5.23p by inversion of the finite Hilbert trans- 
form 



5.3 The planar correlator from the typical tableau 

To deal with the correlator (15.131) . we use the following trick [13]. Fix a tableau R and 
suppose there are a + 1 outward corners. For our purpose there, the useful formula for the 
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Xs = I5 X4 I4 X3 13 X2 Y2 Xi Yi 

Figure 8: Symbols Xj and Yj denote the values of N — j + Rj at the J-th outward and inward 
corners. 



dimension is [TT] 

Ri - Rj + J 



dimi?= n • (5.34) 



l<i<j<N 



If the box [i, Ri) is at the /-the outward corner, the formula fl5.34p gives the residue of the 
pole at = h{N — i + Ri) as 



^ a+l a+1 

j=i j=i 

Here Xj and 1/ denote the value of + Rj — j at the outward and inward corners as shown 
in Fig. [HI Note that Xa+i = 0. We recognize (15.351) as the coefficients in the partial fraction 
expansion of 

a+l 



j=i 

The constant term is fixed by the asymptotics VL{uj) ~ l/cu for large u. Because f"{x) is a 
sum of delta functions at the corners, the second term can be written as 



1 I 

— exp- / dxf" {x)\og{x — UJ IX) 
^ 2 Jo 



Xuj \ 2 Ja_ X — UJ /X 

Here we assumed that /'(x) = — 1 for < x < x_ and that /'(x) = 1 for x > «+. Thus 

n(cj) = - - ±1^ L exp -- / dxf\x) . (5.38) 



A Xuj V 2 x — ujjX) 
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When f{x) is a piecewise linear function representing R, this relation is exact. 

In the large limit, the sum over tableaux is dominated by the saddle point /^,. Thus one 
needs to perform the integral in fl5.38p with / = /* and a± = x± to obtain the correlatorJ§ 
One indeed finds that fl5.38p agrees with the result (12.131) obtained from the SD equation, 
namely, 

G{u) = iQ{uj)\f=f, = '^\~^^ (a + - v/(c^-^+)(^-cu-)) • (5.39) 
Note that u!± = Xx±. 

6 Rectangular models 

The qQ model can be generalized to the case that a]., is a rectangular N x K matrix, while 
the fundamental is still A^-dimensional. For example, in the limit of large N with K/N fixed, 
the graphs of Fig. 2 have weight h'^'^^KN"', and so the SD equation becomes 

• nfrr \ 1 iyhKG{T,uj) 

— iujG(T,uj) = 1 . 6.1 

l-y + ihNG{T,uj) 

The solution is now 



G{T,uj) = ^ [h{N - Ky) + c..(l - y) - {I - yW{uj - uj+){uj - u^) 



h 

y 



{n + yK ±2^NKy^ . (6.2) 



The sum over Young tableaux generalizes as follows, 

-iG{uj) = (1 -i/)^^^^i/l^ldimAri?dimi^i?n(cu) = {n{uj)), (6.3) 

R 

uj-h{N-i + Ri) 

(«,ft,;):corner 

(i,Ri):corncr 

The subscripts A^, K refer to the dimensions of the representations of U{N) and U (K) with 
the given tableaux. The tableaux thus have at most min(A^, K) rows. 

There are several special cases and limits that can be solved more fully. 



For this, change the variable to 6 = {n/2)ip{x), so that the integral takes the form /T^^/j ^ cos6/ (ci 

2 

/2 



sin6')(c2 — sin 6'). It can be integrated by parts using /^^/2 d9\og{c — sin 6') = 7rlog[(c + Vc^ — l)/2] 
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6.1 N = l 



Here, 

H = ha^aAlAi' (6.6) 

is just the product of number operators. The energy is in the zero-fundamental sector and 
hk in the one-fundamental sector, where k is the number of adjoint excitations. Thus, 

oo 
fc=0 

where 

P(fc) = (l-y)V(^ + f (6.8) 
is the normahzed probabihty to find a total of k excitations on K oscillators. Also, 

GW = ^w(r4^)" • (6-9) 

Young tableaux are limited to a single row, with any number k of boxes, so that dim^r R = 1 
and dim^-R = which reproduces the above result. 

Though the exact solution fl6.9p is available, it is interesting to consider a saddle point 
approximation. The dominant tableau at large K is at 

K/{k, + K)=y, k, = Ky/{l-y). (6.10) 

This simply moves the pole of the free propagator from uo = to uj = hKy/{l ~ y). To see 
the continuous spectrum we must look at the fluctuations around the dominant tableau. 



This has a width of order v -ft' so we can replace the sum by an integral, k — k^ xy K 



G{uj)^t -p. ^= . (6.12) 

J-oo v27r7 uj — hKy/ {1 — y) — hyf Kx + ie 

Since this has no singularities in the lower half-plane except the pole, G{t) for t > can 
be obtained by a contour integral in uj. Integrating over x, one finds that the falloff in t is 
gaussian. However, in the exact expression there are periodic recurrences at At = 27r//i. 

We can also look at this graphically, as a large-i^ vector model. One would normally 
hold hK fixed, and the dominant graphs are then tadpoles like the n = term in Fig. 2. In 
our model this is just an uninteresting shift of the frequency, as we see from the result fl6.12p . 
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It is therefore more interesting to absorb this into an additive redefinition of uj and to take 
a different hmit in which C, = h^/K is fixed. The relevant graphs are hke the n = 1 term 
in Fig. 2, where the i^'- vector loop interacts at two points with the fundamental excitation. 
For interaction times t' and t" this loop contributes 

-^(t'-t")CV, (6.13) 

so its integrated contribution is — t^C^7^/2. For = 1 the relative time-ordering of the 
interactions of different loops does not matter, and so the loops exponentiate, 

G{t) = e{t)e-''^'^'/^ . (6.14) 

This is indeed the Fourier transform of the result (16.121) . The falloff is gaussian in the large- 
approximation, while the recurrence time At = 27rV^/C grows as \/K. 

Of course, given the exact solution for = 1 we can write the full 1/K expansion, and 
it it useful to do so. In particular for t > 

InG(t) = i^ln(l -y)- K\n{l - ye*^*/^) 

OO J, 

r=l ^ 



r s\ 

r=l s=l 



For large order s, the dominant terms in the sum are at r ~ s/| ln?/|. Replacing the r sum 
by an integral, the coefficient of K^~'^/'^ is approximately 

{tCt/\\ny\y . (6.16) 



\2 



Thus, the 1/K expansion has a finite radius of convergence, but the radius {m(3/C,t) de- 
creases with t. This is interesting in connection with the recurrences, as we will discuss 
further in the conclusions. 



6.2 Large K, fixed N 

We can extend the graphical solution to give the large-ii' limit for any fixed A^. With 
C, = hy/K fixed, and the tadpole shifted away, the only graphs that survive as Ji" ^ oo are 
bubbles with two vertices, though at finite different bubbles can overlap in a nonplanar 
way. Thus we can replace Qij = A\yA^,- with the U{K) singlet field (pij\fK, with two-point 
function 

(T^,,(t)^,,(t')) = ^ {^Q^,mkl{t')) = -^^M^k . (6.17) 
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Figure 9: a) Loop of A; = 5 vertices and adjoint propagators, b) Dual graph with 
interaction. The arrows point from ip to 0. 



This is independent of time, so ipij is governed by a gaussian matrix (rather than path) 
integral. At fixed ipij the fundamental propagator is -^Tr(e~*^'^*), and so 



For = 1 this is an ordinary integral and reproduces the earlier result. For = 2 we have 
the eigenvalue integral 



For general A^ we can express G{t) using the method of orthogonal polynomials [T6] . 



In general, the adjoint propagators form cycles with k vertices as in Fig. 2, though at finite 
A^ these cycles can overlap, and the time ordering need not match the cyclic ordering. We 
would like to rewrite the cycle in terms a dual graph, as in Fig. 9. The loop becomes a vertex 
Tr-?/;^, attached to the fundamental line by ipip propagators. We must add two fields, with 
an off-diagonal propagator, so that every propagator has one end on the central vertex and 
one on the fermion line. 

In general this seems rather complicated, because we must take account of the number 
of backward and forward propagators (13. ip . However, in the special case y ^ 1, where we 
must hold fixed x = —y), these are equal and the total contribution of the vertices and 




(6.18) 




(6.19) 



6.3 y ^1 limit, all N,K 
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adjoint propagators is simply {—ix)^K. In this case, all graphs can then be summed via 

j(t) = 1J rf^V^^'^Ve-^'-(-*^^+^(^«Tr(e-*^^*) , 
^ (-ill))'' 

V{^) = -KY,^^^ = KHl + t^) . (6.20) 

k=l 

The central vertex contributes K{—i)^, the propagators contribute i'', and the outer vertices 
contribute (-^x)^'- The U{N) group structure of the two graphs is also identical. 

We can check the result for = 1, where it becomes an ordinary integral. The ip integral 
produces 6{ip — xt), and the ip integral then gives e~'^^^'^\ which agrees with the y 1 limit 
of the earlier result (16.91) . Note that at infinite temperature we are accessing an infinite 
number of states, so there will not be recurrences. 

7 Discussion 

The qQ model, though in some ways quite simple, is seen to have a rich expansion. 
Thus it may be a useful laboratory for studying this expansion in connection with AdS/CFT 
duality. In particular, we have presented several ways of analysing the model — graphical, 
loop equations. Young tableaux — each of which is suggestive of a bulk string or geometric 
language. 

Our main interest in this model is that it captures the large- structure of the information 
problem, in that correlators decay to zero at infinite N but not at finite N. It is interesting 
from this point of view to consider the = 1 model of Sec. 6.1, which has the same features: 
there is gaussian decay in the large- limit but recurrences on a time scale \/K. In this 
case we have the whole 1/K expansion, and so we can analyze this further. In fact, from the 
coefficient (16.161) we see that the expansion parameter is proportional to t/^/K. Thus there 
is no surprise that the limits K ^ oo (giving decay) and t ^ (yo (giving recurrences) do not 
commute. 

The existence of recurrences does not require this behavior of the perturbation theory. 
Consider for example the function 

oo 

Gtoyit) = e-(*-^-'v^/^)^^^^^/2 , (7.1) 

/=— oo 

This has the same large- limit as the model, and recurrences on the same time scale, but 
by design there are no perturbative corrections at all. 

The large-i^ limit is a vector model, and there is no reason to expect the same result 
for matrix models. However, it is notable that we have found that the expansion parameter 
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in the matrix model is of order /N"^. We must be careful, however, in interpreting this. 
Consider for example a system having a quasinormal mode e~^*, with V of order T. Quantum 
gravity effects will certainly lead to a small correction to the decay rate, F ^ F + Fi/A^^, so 
that the quasinormal behavior becomes 



Thus, corrections that grow in time (relative to the leading term) are inevitable, even 

for a perfectly well-behaved perturbation expansion]^ Can we distinguish harmless growth 
from growth that leads to recurrences? It may be a clue that the series (16.161) has a finite 
radius of converge, while the series (17.21) converges. Thus it would be extremely interesting 
to determine at least the dominant large-t behavior at each order in to see whether it 

represents a breakdown of the expansion at large times or simply a small shift in the 

leading behavior as above. If there is a long-time breakdown of the l/N"^ expansion in this 
model, it is then important to determine whether it is an artifact of the model, or carries 
over to the full dual gauge theory. 

In the context of real black holes, in order to understand possible corrections it is useful 
to consider Hamiltonian evolution along a series of spacelike surfaces ('nice slices') of small 
curvature (see, e.g., [H]). The absence of any large local curvature invariants then implies 
that gravity as a low energy effective theory does not signal its own break down due to 
high energy effects. The dimensionless loop expansion factor is G/R"^, where Rs is the 
Schwarzschild radius. Dimensional analysis would still allow effects of order Gt/R^, which 
become of order one on the lifetime of the black hole, and such effects were found in [18j. As 
the above discussion shows, it is nontrivial to determine whether such large corrections signal 
a breakdown of the perturbative calculation, or something that is readily resummed. In the 
calculation of Ref. [18] we believe it is the latter: the effect is an incremental change in the 
lifetime of a black hole due to the energy of a Hawking particle, similar to the incremental 
change (17. 2p . Indeed, the no-hair theorem would seem to imply that secular effects can only 
manifest themselves through the mass (and other conserved charges) of the black hole. Thus, 
from the Hawking radiation one can learn something about the amount of mass flowing into 
and out of the black hole over time, but one cannot distinguish two inf ailing objects of the 
same mass but different internal states. 
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^There is nothing particularly gravitational about this example: the same would hold for QED corrections 
to muon decay. 




(7.2) 
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